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1. Introduction 

The finite-temperature plenary talks at the yearly Lattice conferences have traditionally been 
overviews of results obtained during the past year or so [[l]]-[|3|]. This time, using an "outsider" 
status as an excuse, I would like to depart from the tradition and start with a somewhat more 
"active" approach, outlining a few newer ideas with the hope that this may help to inspire future 
work (sec. §). Subsequently, however, I return to time-honoured practices and summarize results 
that caught my attention during recent months (sec. 

More precisely, the basic novelty that I would like to elaborate on is the slight paradigm shift 
that has been taking place concerning the role that heavy quarks (charm and bottom quarks) may 
play in hot QCD (with a temperature T ~ 150 — 500 MeV). Given that the masses of the heavy 
quarks are much above the temperature, by up to an order of magnitude, it was long thought that 
they would be relatively "inert" and play little role in this temperature range. It has been one of 
the remarkable empirical discoveries of the Relativistic Heavy Ion Collider (RHIC) at Brookhaven, 
though, that even within the short lifetime of the thermal system the heavy quarks do appear to 
experience significant interactions with it. In fact, this is one of the reasons why the medium 
generated in heavy ion collisions is nowadays conceived to be a "strongly coupled" one. 

Before embarking on a specific discussion on this topic, let me try to place my presentation 
in a wider context. Indeed, high-temperature QCD, or pure SU(A'^^) gauge theory, can be pursued 
to many different goals. These theories offer, for instance, a tractable theoretical limit in which to 
study various aspects of confinement and chiral symmetry breaking, and several recent papers as 
well as parallel and poster contributions were formulated in this spirit (see, e.g., refs. [Q]). On the 
phenomenological side, the original motivation for considering hot QCD was the possible role that 
it may play in Early Universe cosmology, and indeed the relic density of certain dark matter can- 
didates is sensitive to the QCD equation-of-state (see, e.g., refs. [^). Currently the most pressing 
issue of the field is, however, to offer QCD-based non-perturbative predictions for the observables 
that play a role in on-going and future heavy ion collision experiments, so this will be the focus of 
the present talk. 

2. Heavy quarks at high temperature 

Traditionally it was assumed, based on leading-order weak-coupling computations, that charm 
quarks would not have time to thermalize in heavy ion collisions, neither "kinetically" nor "chemi- 
cally". Kinetic thermalization means that the momenta of the heavy quarks be distributed thermally, 
i.e. that the average momentum vanish in the rest frame of the thermal system and the average 
momentum-squared be proportional to the temperature. Chemical thermalization means that the 
number density (or, more properly, the entropy density) associated with the heavy quarks be as 
large as thermal field theory predicts. 

A kind of a paradigm shift, associated with the concept of a "strongly coupled quark-gluon 
plasma", has however been taking place during the last few years. Indeed, experimental observa- 
tions (to be reviewed below) concerning the "quenching" of heavy quark jets can be interpreted 
as indirect evidence for their kinetic thermalization; that is, heavy quarks interact more strongly 
than originally expected. Whether a chemical equilibration also takes place is not quite as obvious; 
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changes in number density are accompanied by an additional Boltzmann factor, exp(— 2M/r) <^ 1, 
where M is a heavy quark mass and the factor two accounts for the fact that quarks and antiquarks 
come in pairs. On the other hand, the initial state already contains a distribution of heavy quarks 
and antiquarks originating from the very first hard scatterings, and this distribution may happen to 
be of the right order of magnitude; if so, only kinetic thermalization is required in order to bring 
them to full equilibrium. In any case it is now a challenge for theorists to, first of all, understand 
quantitatively the rapid kinetic thermalization of heavy quarks, and second, just in case, to probe 
how big an effect chemically thermalized charm quarks would have in the hydrodynamic modelling 
of heavy ion collision experiments. 

In order to organize the corresponding discussion, I start by considering the case of "0 valence" 
heavy quarks; by this I refer to the effect of (chemically thermalized) heavy sea quarks on the 
equation of state. This case can also be identified as QCD with Nf = 2 + \ + \ flavours. I then 
proceed to the "1 valence" sector, considering heavy quark jets and their kinetic thermalization; 
and end with "1 + I valence" heavy quarks, meaning heavy quarkonium. As we will see, there is a 
perspective for progress on all these fronts. 

Before proceeding, I would like to briefly ponder the question of when it is precisely that some 
quark is "heavy", compared with the temperature T . To this effect, recall the form of the free quark 
propagator in continuum: 



with (On = ±7iT, ±37iT, .... We see that the relevant comparison is something like M ^ kT , but a 
question remains whether, once interactions are taken into account, we should insert the MS mass, 
say Mf\3 GeV) « 1 GeV; the pole mass, say Mr'*" ~ (1.5 - 2.0) GeV; or something else. There 
is probably no unique answer to this question; rather, the answer depends on the observable. In 
any case, the message to take home is that in principle charm quarks could be "light" as soon as 
T > 1 GeV/TT ~ 300 MeV, or "heavy" as long as J < 2 GeV/Ti ~ 600 MeV. 

2.1 "0 valence" — charm quark effect on {e — 'ip)/T^ 

It was suggested a few years ago, based on a next-to-leading order weak-coupling analysis, that 
charm quarks may have a significant effect on various thermodynamic observables at surprisingly 
low temperatures [^]. The issue then is whether lattice studies could consolidate this suggestion. 

As far as the lattice goes, it is appropriate to point out that although we conceptually now 
consider QCD with A^f = 2 + 1 + 1 flavours, with charm quarks in the thermal sea, the lattice 
estimate of their effect can in practice be reduced to the measurement of certain condensates, which 
can then be evaluated within the 2 + 1 flavour theory. In this approach the charm quarks are only 
partly dynamical. Nevertheless, the order of magnitude of their effect should still come out right. 

More concretely, focussing on the trace of the energy-momentum tensor, i.e. {e — 3p) /T^, the 
heavy quark contribution can be reduced to |0] 



where Zi oc dM^a^e/dpi^ is a coefficient function that needs to be determined along the lines of con- 
stant physics (/3l refers to the lattice j3 -coefficient). Recent measurements, from asqtad fermions 
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Figure 1: The charm quark contribution to the trace anomaly according to recent lattice simulations (left) | 
as well as resummed perturbation theory (right) [g] . 



withA^T = 6 and bare mass ratios Mbare/'«i ~ m^rriu ^ 10 [||], are compared with the corresponding 
curve from perturbation theory in fig. |l|. 

We observe that while the effect of charm quarks is very small at T < 250 MeV, their relative 
importance in the trace anomaly increases rapidly with temperature, reaching >50% already at 
r>400 MeV. (The relative effect is somewhat smaller in quantities like the pressure or energy 
density, where no subtraction is carried out.) 

To summarize, if chemical equilibration takes place, then charm quarks might affect the initial 
stages of hydrodynamics in future heavy ion collisions at the LHC, where higher temperatures 
may be reached than at the RHIC. In any case, the charm quarks do play a significant role in the 
equation of state relevant for cosmology, in which enviroiunent chemical equilibrium is guaranteed 
to be reached for all strongly interacting particles. 



2.2 "1 valence" — heavy quark jets 

In any high energy collision, a number of heavy quarks (and antiquarks) are produced in an 
initial hard process, as can be illustrated by the following Feynman diagram: 




Assuming a suitable factorized framework and making use of previous experience gathered in p+p 
or p + p collisions, it is believed that the corresponding production cross section is relatively well 
understood also for collisions involving heavy nuclei, such as (i + Am and Am + Am Subsequent 
to their production, the heavy quarks decay, often semi-leptonically as c ^ IvX. The leptons £ can 
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Figure 2: Experimental results from ref. [ p^ (see also ref. [|ri]|), indicating that heavy quark jets propagate 
as expected ind + Au collisions, but get quenched in Au + Au collisions. 



be observed, and the outcome can then be compared with the theoretical prediction. 

It turns out that while results from d + Au collisions indeed conform with theoretical expec- 
tations, those from Am + Au collisions appear not to do so; rather, less leptons i are observed than 
expected. The results are illustrated in fig. Bin terms of the so-called "nuclear modification factor". 



Raa, which is significantly below unity for Am -i- Am collisions [|10p; we can say that heavy quark jets 
get "quenched". Another relevant observable is the so-called elliptic flow, and the indication is that 
heavy quarks do participate in this hydrodynamic behaviour [|ri|]. Both observations point towards 
the interpretation that, due to multiple scatterings with other particles, the heavy quarks slow down 
with respect to the thermal medium, and then flow together with it. I refer to this phenomenon 
as kinetic thermalization. In fact, the heavy quarks behave much like heavy particles in classical 



non-relativistic Brownian motion, and many of the same concepts can be argued to apply []12|]. 

Let us denote the rate at which heavy quarks, assumed already to be close to rest, slow down, 
by r]D; this thermalization rate is often also called the "drag coefficient". Inspired by the anal- 
ogy with Brownian motion and its classical description through Langevin dynamics, r\j) can be 
fluctuation-dissipation-related to another coefficient, fc, characterizing the autocorrelation of the 



force that acts on the heavy quarks [|12|, [13|] : 



JC / a^J^rW ,. 2Tpe((o) 

riD = 1 + — , k: = hm — !l£l—L (23) 

2Mkinr \ \Mi,in J J cn^O (0 ' 

where M^in refers to a particular heavy quark mass definition, related to the pole mass at zero 
temperature. The function p£ is the spectral function corresponding to the Euclidean correlator 
(assuming T > or Nf > 0) ^ 



G£(T) = --£ 



I 3 (ReTr 



Up.,,gEiiT,0)U,.flgEi{0,0) 



i=i {ReTr[Up.,o 



(2.4) 
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Figure 3: Results for the force-force transport coefficient K within classical lattice gauge theory, plotted as a 
function of the perturbative expansion parameter g^N^T /niD i^^i (from ref. [|l4|]). The non-perturbative data 
can exceed the leading-order and even the next-to-leading order results by a significant amount. 



where gEi = j[Do,D,] is the colour-electric field (shown here in continuum notation), and Uz^^,z.^ is 
a Wilson line in the Euclidean time direction, from to Tb. The numerator can be illustrated as 




where the circle represents the Polyakov loop around the Euclidean time direction and the blobs 
denote electric field insertions. As far as I know no lattice measurements of this correlator have 
been published yet, although they should not be overwhelmingly demanding. 

Even though no actual measurements are available, another test has been carried out. Indeed, 
the electric field correlator can also be addressed within "classical lattice gauge theory", which 
serves as kind of an effective low-energy description of infrared phenomena at finite tempera- 



tures [|15p. It turns out that classical lattice gauge theory is more sensitive to ultraviolet physics 
than would be expected from a proper effective field theory framework; nevertheless, the result is 
interesting if plotted in terms of a quantity that does have a direct analogue in QCD, the so-called 
Debye mass parameter, m^^^^^ ~ g^T /a. A result is shown in fig. ^ including comparisons with a 
leading order perturbative result and a fitted next-to-leading order behaviour (the next-to-leading 
order result has been computed analytically in continuum QCD [|^], but not in classical lattice 
gauge theory, where only its parametric form is known). 

It can now be observed that if we insert the estimate g^T ~ n?D,iatt> which is known to be a 
reasonable one on the QCD side, then we are in a regime where the next-to-leading order correction 
is of order 100%, yet the non-perturbative result is larger still. Excitingly, such a significant increase 
appears to be more or less on top of what is phenomenologically needed in order to explain the fast 
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Figure 4: The combination 2T Pe{(o) / (O from eq. (2.3), in lattice units, according to classical lattice gauge 
theory (from ref. [[l4|]). The intercept at O = yields the result plotted in fig. ^ in different units. 



thermalization (quenching) of heavy quark jets (see, e.g., ref. [|17|] and references therein). 

Classical lattice gauge theory makes another prediction as well. In fig. ^ the frequency- 
dependent function from which the intercept is to be taken according to eq. (^j|), is shown. The 
basic observation is that this function is/?af at small frequencies; it has no transport peak, unlike 
spectral functions related to conserved currents. A similar flat behaviour for the spectral function 
related to the electric field correlator has also been observed in a very different theory, strongly cou- 
pled ,jy = A Super- Yang-Mills theory in the large-A'^c limit, handled through its AdS/CFT dual [^]. 

All this suggests that the electric field correlator might be more amenable to analytic continua- 
tion than current-current correlators from which viscosities and conductivities have been extracted 
previously [p^. To summarize, I can only solicit numerical tests of eq. (2.4), and hope that the 
outcome might yield a large coefficient K, differing from the perturbative one by as much as an 
order of magnitude. With the advent of LHC, which should produce data on bottom quark jets as 
well, containing a different M^in in eq. (13) than in the charm quark case, our understanding of 
heavy quark jets within a hot medium could then be quantitatively tested. 



2.3 "1 + 1 valence" — heavy quarkonium 

As a last example of heavy quark related observables, I briefly summarize recent news from 
heavy quarkonium physics. Like heavy quarks, heavy quarkonium can originally be generated in a 
non-thermal hard scattering, or through a thermal fluctuation. The latter process is reminiscent of 
those contributing to the chemical thermalization of single heavy quarks, and therefore probably 
too slow to take place effectively; nevertheless, surprises cannot be excluded. After having formed 
one way or the other, heavy quarkonium propagates through the thermal medium, whereby its prop- 
erties get modified; therefore the quarkonium peak observed in the dilepton rate [^] may change 
in magnitude, shape, or position, depending on the temperature that is reached in the collision [pT|]. 
Some relevant Feynman diagrams are illustrated in fig. ^. 
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Figure 5: Left: an amplitude corresponding to the production of a dilepton pair from a thermalized medium. 
Right: squaring the ampHtude, the dilepton production rate is seen to be proportional to the two-point cor- 
relation function of the electromagnetic current, with significant near-threshold corrections coming from 
medium-modified Coulomb exchange. 



Now, like always in quarkonium physics, it appears reasonable to try and address the thermal 
modification of quarkonium properties through a potential model, thereby resumming corrections 
from graphs of the type in fig. fright). At finite temperatures, however, this is complicated by 
the multitude of different potentials that can in principle be defined, and intuitive arguments alone 
cannot decide which of them is the correct one. What is needed is rather a derivation of the relevant 
effective framework from QCD; at least within perturbation theory this can indeed be achieved [^], 
resulting in a potential-model type recipe for computing the spectral function [p3|], with a definite 
(in general complex) potential appearing as a "matching coefficient". 

In order to appreciate the intricacies of the issue, it is important to realize that, on a Euclidean 
lattice, the time extent j8 = 1/r is in some sense always "small"; more quantitatively, it can be ar- 
gued that quarkonium melts at a temperature where parametrically jS < l/a^Mkin p3|]. In contrast, 
the Minkowskian time scale t corresponding to the Coulombian binding energy of heavy quarko- 
nium is "large", f ~ 1/ a^M^^ai- So, in the heavy quark hmit where the effective is small, we see 
that f S> j8 , and it is more or less clear that the potential relevant for discussing quarkonium binding 
and dissociation at finite temperatures should involve some sort of an analytic continuation. 

In this conference, a very interesting suggestion for a non-perturbative definition of a real- 



time static potential in this spirit was put forward [|24j]. Motivated by the perturbative definition, 
the idea is to first measure a Wilson loop as a function of a Euclidean time coordinate, t; this can 
be illustrated as 



which observable we denote by C£(T,r) = (Tr[W£(T,r)]). Supposing that an analytic continuation 
can be carried out, T it, a real-time potential could then be extracted from 

id,CE{it,r) = V>{t,r)CE{it,r) . 

Its static limit corresponds to y>(oo,r). 

It is appropriate to point out that the horizontal Wilson lines appearing in the definition of 
C£'(T,r) are non-unique, as usual; however, to the extent that we can compensate for the specific 

1/2 ~ 

choice by a normalization factor, (r), so that the correlator reads C£(T,r) = Zv(r)C£:(T,r), 
we see that Zv(r) drops out from the definition of Vy{t,r). Such a normalization factor is indeed 
characteristic of the effective field theory framework that can be used for addressing the properties 
of heavy quarkonium 
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Figure 6: Left: the trace anomaly (e — 3p)/T^ fromref. [p6|], with two famihes of fermion discretizations. 
Judging with bare eye, the results appear to converge towards a common continuum limit. Right: baryon 
number, electric charge, and strangeness susceptibilities (Xi' Xi^ respectively), from ref. [p7|]. 



Now, according to ref. [24], the static limit y>(oo,r) can indeed be extracted through a spec- 
tral analysis of the Euclidean correlator. More precisely, the data appear to indicate the pres- 
ence of a spectral peak, whose position signals the average energy of the quark-antiquark system, 
ReV>(oo,r). The peak should have a finite width as well, ImV>(oo,r), being a signal of a Coulomb 



scattering/Landau damping induced "decoherence" of the quark-antiquark state []22|], caused by 
collisions with the particles of the thermal medium. 

To me, the idea of ref. [^J] seems very interesting, and I am looking forward to further devel- 
opments along these lines. 



3. Light quarks and gluons at high temperature 

I now move away from heavy quark related observables and discuss a number of recent devel- 
opments related to light quarks and gluons. 

3.1 Basic thermodynamics witli Nf = 2 + \,Nc = 3 

Results for the trace of the energy-momentum tensor (like in fig. |T]) and for various physical 
susceptibilities, from the large-scale simulations in refs. [ p6| , p7| ], are shown in fig. ^ These re- 
sults are supposed to be "physical", i.e. for (almost) realistic quark masses, and the lattice spacing 
dependence would appear to be regular as well. Improved staggered quarks were used. 

It is fair to say, though, that it is not easy to judge the systematic uncertainties that may still 
be hidden in these results. For instance, in fig. ^ a comparison of the strangeness susceptibilities, 
as determined by the RBC-Bielefeld and by the Budapest-Wuppertal collaborations, is shown [^]. 
The infamous 30 MeV temperature shift is clearly visible. 

The good news is that the observable of fig. ^is a nice one: it is directly physical, being related 
to the zero component of a conserved current, so that no renormalization ambiguities should appear. 
I encourage other groups to also primarily carry out comparisons with this quantity, rather than with 
more ultraviolet sensitive "auxiliary" observables, like the susceptibilities related to the Polyakov 
loop or the chiral condensate. 
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Figure 7: The strangeness susceptibilities from different collaborations, according to ref. 



As far as the reason for the discrepancy goes, I am certainly not qualified to offer any original 



insight. I do assume that "trivial" issues, like improper thermalization (see, e.g., ref. [pOl]), have 
been carefully excluded by all collaborations. As a "user" I do find it a pity, though, that none of 
the collaborations shows finite- volume scaling in their plots, even though the chiral limit is not far 
and the transition is very weakly of the first order; indeed figs. ^, ^ refer to a fixed box size L in 
units of the temperature, LT = 3 — 4. ' On the side of discretization effects, it would be nice to 
overlay results from Wilson-like discretizations on fig. ^, particularly given the delicate role the 
chiral symmetry plays close to the transition point, and fortunately efforts in this direction appear 



to be under way pi|, 32, 33]. 



3.2 New precision for Nf = 0, Nc = 3 

Recently, there have been new studies of the thermodynamics of pure SU(3) gauge theory. The 
great benefit of this simplified theory is that systematic errors can be brought better under control; 
therefore the theory offers an excellent test bench both for new lattice ideas, and for comparing the 
lattice data with various continuum computations. These two aspects are illustrated in fig. |[ On the 



left, a test is shown of the new approach of ref. [|34|], in which the lattice spacing is kept fixed and 
the temperature is varied through changing the number of points in the Euclidean time direction, 
N^; this is theoretically more transparent than the standard approach where A^^ is kept fixed and 
temperature is varied through j3L, implying a simultaneous variation of the lattice spacing. On the 
right, a very precise study of the entropy density at low temperatures is shown; the entropy density 
is a convenient observable in that it can be measured without any subtractions (because the entropy 
density of the vacuum state vanishes), through 



3T 



Z2Tr[B2-E2 



(3.1) 



where is a renormalization factor. It can be seen that the results are precise enough to allow for 
a stringent comparison with the contributions from various glueball spectra. 



Very recently scaling studies in the critical region have been reported in ref. | 
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Figure 8: Left: the trace anomaly (e — 3p)/T^ from ref. [ p4j | for A^f = 0, in the so-called fixed-scale ap- 
proach. Right: the entropy density from ref. [^Sjl for A'f — 0, together with a comparison with various glueball 
resonance descriptions. 



In my opinion, this kind of precision relooks at pure gauge theory are very welcome, and 
deserve to be pursued for many other observables as well. 

3.3 Another look at A^^f = 0, A^^c > 3 

Building on the previous section, it is also interesting to inspect pure gauge theories with A^c > 
3. Recently two groups have come up with new results in this spirit (cf. fig. In ref. [|3^], results 
have been presented for A'^c = 3,4,6, with a focus on scale setting and continuum extrapolation 
{Nx = 6, 8). Another work was presented in ref. [^] where, following the earlier work in ref. [Q], 
results were presented for A'c = 3, 4, 5, 6, 8, for a fixed = 5. 

A very intriguing trend can be extracted from the results of ref. [ |3^ ] (fig. ^(right)). Indeed, 
it appears that for large A'^^, the functional form of the trace anomaly becomes much simpler than 
for A^c = 3: once normalized to N^T^, the result is basically zero for T kT^; displays a large jump 
(a first order transition) at T = Tc; and decreases then monotonically for T > T^. This is a very 
simple pattern which suggests, and even calls for, a theoretical explanation; in fact, it can perhaps 
be speculated that in this limit the high-temperature phase is qualitatively purely perturbative. 

A word of caution may be in order, though. In ref. [09|], the latent heat (the discontinuity of 
the energy density or, for that matter, of the trace anomaly) was studied at various A'^^ for A^^ = 
5,6,8, and it was found that in general the results at A^^ = 5 did not fit well into the continuum 



extrapolation, particularly at large N^- Given that in ref. [ p7[ ] only Nt = 5 was considered, the results 
should probably be assigned a generous systematic error (as indeed elaborated upon in ref. [^); 
in particular, it might be desirable to have additional A^^'s to add confidence to the Nc ^ °° limit 



(for A'^c = 4, 6, the results of ref. [36] might be helpful in this exercise). In any case, 1 consider this 



to be a very interesting topic and worth further study. 



3.4 Finite-volume effects and screening masses 



In section p| , 1 already alluded to finite-volume effects. Recently, an interesting theoreti- 
cal work appeared [ pO[ ] in which the finite-volume effects in various thermodynamic observables 
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were analyzed in some detail. For instance, for the entropy density, s = {p-\-e)/T , which can be 



measured according to eq. ( |3.1| ), the expression 



-m{T)L o , . 

s{T,L) = s{T,oo) - V-— — [1 + ^TdT]m\T) + ^(e-v^^OT^) (3.2) 



2kL ' 2 



was given, with m{T) denoting the lightest screening mass and v its degeneracy. Similar formulae 
were also given for combinations like e — 3p and p; the ones for e — 3/? would play a role in mea- 
surements such as those in figs, ^left), ||(left), ^(right). The exponential dependence in eq. (3.2) 



is familiar, but it is interesting that the pre-exponential factor is also completely fixed in terms of 
known or measurable quantities. 

In my opinion, the fact that finite-volume effects are determined by the lightest screening 
mass, as exemplified by eq. ( [3.2| ), underlines the usefulness of measuring the screening masses in 
the context of every lattice study. Indeed, though of no direct use for heavy ion experimentalists, 
the screening masses allow theorists to learn a lot about the dynamics of the system, and also to 
judge whether systematic errors related to finite- volume effects can be under control. 

Examples of measurements of screening masses are shown in fig. |l^. Among the qualitative 
observations that can be made from these plots are the following: 

• For T > l.5Tc, the screening masses coupling dominantly to gluonic objects (i.e. flavour sin- 
glets) are the lightest ones. That this should happen is one of the predictions of the dimen- 



sionally reduced effective description of high-temperature QCD [|43|]; thus, we may assume 
that dimensional reduction could work qualitatively for T>l.5Tc while it probably cannot 
capture all the relevant dynamics for T < l.STc if physical Nc and A'f are being used. 

It can be seen that cutoff effects in the mesonic screening mass become substantial at high 
temperature. It is interesting to compare this behaviour with that in fig. 0; much the same 
appears to be the case in the strangeness susceptibility. 

At low temperatures, on the other hand, cutoff effects in the mesonic screening mass appear 
to be moderate, while at the same time the screening mass itself becomes small, m/r ~ 
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Figure 10: Left: the smallest screening masses extracted from glueball like operators, with Nf — 2 Wilson 
quarks in a V = 16^^ x 4 volume (from ref. pl|]). Right: the smallest screening mass extracted from a 
pseudoscalar mesonic operator, with A^f = 2 + 1 p4fat3 staggered quarks (from ref. [H)- 



1 around T ~ Tc, indicating the "vicinity" of a second order transition. This suggests, in 



accordance with the discussion around eq. ( |3.2D , th&t finite-volume ejfects may be significant 
at low temperatures. For pure SU(3), this is clearly visible also in fig. §(left). 

Hopefully these points serve to illustrate that there is much to learn from screening masses, and 
encourage them to be adopted as a standard part of every finite temperature lattice study. 

3.5 Energy-momentum correlators and sum rules 

I wish to end with a somewhat amusing recent episode, related to the determination of transport 



coefficients. As was already discussed in section |2.2| , extracting a transport coefficient from a 
conserved current is very difficult in general, because there may be a narrow transport peak in the 
corresponding spectral function, p{(o), around ft) = 0, and one would need to extract the height of 
the peak. If, however, analytic information is available on the shape of the spectral function, then 
sum rules relate integrals over p to various thermodynamic quantities, which are easier to measure. 

Now, a particular quantity considered in this spirit is the standard hydrodynamic transport 
coefficient known as bulk viscosity. The sum rule states that 



(0 



P?l'o(«) 



(3.3) 



where 5 is a certain local thermodynamic observable. Setting aside the question (which I consider 
to be very difficult) of whether the functional form of p' 



.bulk/ 



(0) 



Pr = o(^) ^^"^ be reasonably 
modelled, there has even been confusion about how the left-hand side, S, looks like. In fact, two 



different answers have been proposed: one in refs. ||44|, ^], and more recently another one in 
ref. [ p^ ] (see ref. [47] for a related discussion). The two expressions are compared in fig. 11. 



The reason for the difference is in some sense subtle, though eventually perhaps also simple to 



understand []46[]. The older expression applies in Euclidean spacetime; it gives the value of a certain 
two-point function at a fixed Matsubara frequency ft),, = 0, extrapolated subsequently to vanishing 
spatial momentum. The other expression, in contrast, assumes first the limit of vanishing spatial 
momentum, while still keeping full information on all Matsubara frequencies ft),,, in order to allow 
for an analytic continuation to Minkowskian frequency CO. 
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Figure 11: A comparison of two different quantities that have been thought to be relevant for the bulk 
viscosity sum rule, eq. (3.3) (from ref. [Sq]). The lattice data is for pure SU(3) and originates from ref. pq]. 



Clearly, fig. [11| demonstrates that it is important to have the correct left-hand side in the sum 



rule of eq. ( |3.3| ), if information concerning the right-hand side is to be extracted. 



4. Conclusions 



The main point which I wanted to illustrate in this talk is that, apart from the obvious need 
to continue numerical efforts in order to reach chiral, infinite-volume, and continuum limits with 
controlled systematic errors for physical QCD, there is also room and even need for various types 
of theoretical contributions to the understanding of strong interactions at high temperatures. This 



includes both analytic computations, examples of which were given in sections 2.2, 3.4 and 3.5 



as well as numerical efforts in simplified theories, in which the systematic errors can be controlled 



more easily than in full QCD; examples were mentioned in sections [2.3| , p.2| , |3.3[ Hopefully more 
groups will join also these physically motivated "low-cost" efforts in the coming years! 
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